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^~^ • Abstract: In the minimal formulation of gravity with Lifshitz- type anisotropic scaling, 

^\j ■ the gauge symmetries of the system are foliation-preserving diffeomorphisms of spacetime. 

j~^ . Consequently, compared to general relativity, the spectrum contains an extra scalar graviton 

^? , polarization. Here we investigate the possibility of extending the gauge group by a local U{1) 

symmetry to "nonrelativistic general covariance." This extended gauge symmetry eliminates 
the scalar graviton, and forces the coupling constant A in the kinetic term of the minimal 

k^ ^ formulation to take its relativistic value, A = 1. The resulting theory exhibits anisotropic 

^ ■ scaling at short distances, and reproduces many features of general relativity at long distances. 
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1. Introduction 

The idea of gravity with anisotropic scaling [1-3] has attracted a lot of attention recently. 
There are two, somewhat distinct, motivations for developing this approach to gravity. The 
first is driven by the long-standing search for a theoretical framework in which the classical 



theory of gravity is reconciled with the laws of quantum mechanics. A successful outcome of 
this search would result in a mathematically self-consistent framework for quantum gravity, 
not necessarily subjected to experimental tests. Examples already exist - the ten-dimensional 
super symmetric vacua of string theory belong to this category. The second motivation comes 
from a goal which is more narrow, and also much more ambitious: To find, within such a self- 
consistent quantum gravity framework, a theory that reproduces the observed gravitational 
phenomena in our universe of 3 + 1 macroscopic dimensions. 

Both of these motivations are relevant for the development of gravity with anisotropic 
scaling. For a large class of possible applications, it does not matter whether or not the theory 
matches general relativity at long distances, or conforms to the available experimental tests 
of gravity in 3 + 1 dimensions. A mathematically consistent quantum gravity which lacks this 
phenomenological matching can still be useful in the context of AdS/CFT correspondence, 
and produce novel gravity duals for a broader class of field theories, of interest for example in 
condensed matter applications. It can also have interesting mathematical implications, given 
the close connection between the theory formulated in [1,2] and the mathematical theory of 
the Ricci flow on Riemannian manifolds. 

However, explaining the observed features of gravity in our universe of 3 + 1 dimensions 
is still perhaps the leading motivation for developing a quantum theory of gravity. Therefore, 
it makes sense to ask how close we can get, in the new framework of gravity with anisotropic 
scaling, to reproducing general relativity in the range of scales where the laws of gravity have 
been experimentally tested. 

The comparison to general relativity is facilitated by the fact that in the framework 
proposed in [1-3], gravity is also described simply as a field theory of the dynamical metric 
on spacetime. Unlike in general relativity, however, the spacetime manifold M (which we take 
to be of a general dimension D + 1) is equipped with a preferred structure of a codimension- 
one foliation T by slices of constant time, S(t).^ In the minimal realization of the theory, 
reviewed in Section ^]^, the gauge symmetries of the system are the foliation-preserving 
diffeomorphisms Diff(M, J^). Since this symmetry contains one less gauge invariance per 
spacetime point compared to the full spacetime diffeomorphisms Diff(M), the spectrum of 
the linearized theory around flat spacetime contains one additional, scalar polarization of the 
graviton. 

At short distances, the anisotropy between space and time is measured by a nontrivial 
dynamical critical exponent z > 1, leading to an improved ultraviolet behavior of the theory. 
At long distances, on the other hand, the theory is driven to an infrared regime where it shares 
many features with general relativity. First of all, under the influence of relevant terms in the 
classical action, the scaling becomes naturally isotropic, with the relativistic value of z = 1. 
Moreover, the lowest-dimension terms that dominate the action in the infrared are exactly 
those that appear in the ADM decomposition [4] of the Einstein-Hilbert action: The scalar 



^For simplicity, we will assume throughout this paper that the leaves E(f) of this foliation all have the same 
topology of a D-dimensional manifold E. 



curvature term, which sets the value of the effective Newton constant, and the cosmological 
constant term. 

Thus, in the low-energy regime, the action of the minimal theory with anisotropic scaling 
looks very similar to that of general relativity. However, this similarity has its limits, and the 
theories are clearly different even in the infrared. The differences can be understood in three 
related ways: As a difference in gauge symmetries, a difference in the graviton spectrum, 
and a difference in the number of independent coupling constants. First, in the theory with 
anisotropic scaling, the gauge symmetry is reduced to Diff (M, J-"), and the theory propagates 
an extra scalar polarization of the graviton. In addition, the kinetic term in the action allows 
for an additional coupling A, which is undetermined by any symmetry of the minimal theory, 
and therefore expected to run with the scale in the quantum theory. In general relativity, 
the spacetime diffeomorphism symmetries force A = 1 and protect this value from quantum 
corrections. Stringent experimental limits on the value of A have been advocated in the 
literature [5,6], suggesting that at least in this class of models, it must be very near its 
relativistic value A = 1. However, in the regime near A = 1, difficulties with the dynamics of 
the additional scalar graviton have been pointed out (see, for example, [7-9]). 

In order to get closer to general relativity, it is tempting to focus on the structure of 
gauge symmetries. However, this needs to be done cautiously, keeping in mind that gauge 
symmetries are just convenient redundancies in the description of a physical system, and 
therefore to some extent in the eye of the beholder. The more physical perspective is to focus 
on the spectrum of propagating degrees of freedom.^ Thus, we will be interested in finding an 
extension of the gauge symmetry that will turn the extra, scalar polarization of the graviton 
into a gauge artifact. A second option would be to find a mechanism for generating a finite 
mass gap for the scalar graviton - in this paper, we concentrate on the first possibility. 

We will find such an extended gauge symmetry, with as many generators per spacetime 
point as in general relativity. This gauge symmetry can be viewed as representing "non- 
relativistic general covariance" in gravity with anisotropic scaling. The extended symmetry 
eliminates the scalar polarization of the graviton from the spectrum. As a bonus, we find that 
the extended gauge symmetry requires A = 1, thereby reducing the kinetic term to coincide 
with that of general relativity. It is in fact important that our entire construction depends 
only on the form of the kinetic term, and therefore does not restrict the form of the potential 
term in the action. Hence, at short distances, the covariant theory can exhibit the same 
improved ultraviolet behavior associated with z > 1 in the minimal theory of [1,2]. 



^Of course, one can make the theory diffeomorphism invariant in a trivial way, by integrating in the gauge 
invariance without changing the number of physical degrees of freedom. This leads to a theory which is 
formally generally covariant, but effectively equivalent to the original model. Perhaps any sensible theory 
can be "parametrized" in this way [10], and formally rewritten as a generally covariant theory. This process 
of covariantizing a given theory is very closely connected to the Stiickelberg mechanism used prominently in 
particle physics (see, e.g., [11] for a review). It has been applied to the models of gravity with anisotropic 
scaling [1,2] in [7,12,13]. In this paper, instead, we are interested in the nontrivial extension of gauge symmetry 
which actually reduces the number of physical degrees of freedom. Investigating how the resulting theory then 
responds to the Stiickelberg trick is an interesting question, beyond the scope of this paper. 



In this paper, our perspective is that of (effective) quantum field theory, but we restrict 
our analysis to the leading tree-level, or classical, approximation. Quantum corrections are 
expected to modify the scaling behavior of our models, but they are beyond the scope of this 
paper. In addition, our analysis will be strictly local: We freely integrate by parts and ignore 
total derivative terms. Boundary terms play a notoriously central role in relativistic gravity; 
it will be important to extend our analysis to include their precise structure and clarify their 
role in theories of gravity with anisotropic scaling. These are among the many interesting 
issues left for future work. 

The main result of the paper is the construction of the generally covariant gravity with 
anisotropic scaling, which we present in Section ^. Sections |l], ^ and ^ prepare the ground for 
a better understanding of the main results, and explore a few additional issues of interest. 

1.1 General covariance 

In order to explain what exactly we mean by "general covariance," we first consider two 
theories - general relativity, and the ultralocal theory of gravity [14, 15] - and illustrate our 
point using the Hamiltonian formulation of these two theories. 

In fact, throughout this paper we will often resort to the Hamiltonian formalism,^ for a 
number of reasons. First of all, the time versus space split of the Hamiltonian formalism is 
particularly natural for gravity with anisotropic scaling. More importantly, the technology 
available in the Hamiltonian formalism allows us to get a precise count of the number of 
propagating degrees of freedom, and offers a better insight into the structure of the gauge 
symmetries of the theory. Indeed, one of the advantages of the Hamiltonian formulation 
is that one does not have to specify the gauge symmetries a priori. Instead, the structure 
of the Hamiltonian constraints provides an essentially algorithmic way in which the correct 
gauge symmetry structure is determined automatically [10]. In the process, the consistency 
of the equations of motion is tied to the closure of the constraint algebra and the preservation 
of the constraints under the time evolution. Once the full system of constraints has been 
determined, the constraints are separated into first-class (whose commutators with other 
constraints vanish on the constraint surface) and second-class (whose commutators define a 
nondegenerate symplectic form). As an additional benefit, after determining the numbers Ci 
of first-class and C2 of second-class constraints, the number of degrees of freedom N can be 
reliably evaluated by the standard formula [10] 

AA=^(dimP-2Ci-C2), (1.1) 

where d\mV is the number of fields in the canonical formulation [i.e., the dimension of phase 
space). In local field theory, this formula can be interpreted per spacetime point, giving 
the number of local degrees of freedom. We will use this formula repeatedly throughout the 
paper. 



''For the canonical reference on Hamiltonian systems with constraints, see [10]. 



In canonical general relativity [4, 16, 17] on a spacetime manifold M with D + 1 coor- 
dinates (x*,t), the algebra of constraints contains the "superhamiltonian" ?^_l(x) and the 
"supermomentum" '?^j(x), and the total Hamiltonian is just a sum of constraints: 

H= f d^x [NU^ + N'Ui) . (1.2) 

N and N"^ are the lapse and shift variables of the metric, and T-L^^ and T-Li are functions of the 
spatial components gij of the metric and their canonically conjugate momenta vr*-'. Since the 
constraints are all first-class, they generate gauge symmetries, whose generators are 



nc) = y d^xr(x,t)i^,(x,i), n^ie) = y d^xcO(x,t)i^^(x,t). 



(1.3) 



Their commutation relations are well-understood, even though they do not quite yield the 
naively expected spacetime diffeomorphism algebra. True, the commutator of two T-LiS 

me)Mc')] = neokC - c'^^r) (1.4) 

reproduces the algebra of spatial diffeomorphisms Diff(S), and 

[^(r),^±(c°)] = ^±(e'5fcC°) (1.5) 

just states that 'H± transforms correctly under Diff(S). However, the commutator of 'H± 
with itself gives a field-dependent result, 

[^±(c°),^±(c°)] = -<Tn {g^Kedj^ - ed,e)) . (i.e) 

Here a denotes the signature of spacetime: a = —1 for general relativity in Minkowski 
signature. This generalized, Dirac algebra is the Hamiltonian manifestation of the original 
diffeomorphism symmetry (and general covariance) of general relativity, with D + 1 gauge 
symmetries per spacetime point. 

Our second example is the ultralocal theory of gravity, which results from dropping the 
spatial scalar curvature term R in the action of general relativity. Of course, this step selects 
a preferred foliation T of spacetime, and therefore violates spacetime diffeomorphism invari- 
ance. One might naively assume that the symmetry is reduced to the foliation-preserving 
diffeomorphisms Diff(M, J^). However, the analysis of Hamiltonian constraints reveals a sur- 
prising fact [14, 15]: The theory is still gauge invariant under as many gauge symmetries per 
spacetime point as general relativity. In contrast with general relativity, the ultralocal theory 
exhibits a contracted version of the Hamiltonian constraint algebra, with (^]^) replaced by 
its o" —7- limit: 

[n^{i\ne)] = ^, (1.7) 

while the remaining commutators ( |1.4p and ( |1.5D stay the same. 

The theory is "generally covariant" - it has the same number D + 1 oi (nontrivial) 
local gauge symmetries as general relativity, even though the algebra in the o" — )• limit 



still preserves the preferred spacetime foliation structure. Effectively, the spatial diffeomor- 
phism symmetries have been kept intact, but the time reparametrization symmetry has been 
linearized, and its algebra contracted to a local U{1) gauge symmetry.^ Just as the Dirac al- 
gebra of Hamiltonian constraints ( |1.4| ), (1.5) and ( |1.6D in general relativity is associated with 



the Lagrangian symmetries described by the group of spacetime diffeomorphisms, Diff(M), 
the Teitelboim-Henneaux algebra (|1.4D, (1.5) and (p^) can be associated with a Lagrangian 



symmetry group which takes the form of a semi-direct product, 

?7(1) KDiff(Af,J^). (1.8) 



It is natural to interpret ( |l.8| ) as the symmetry group of "nonrelativistic general covariance." 
This is the structure of gauge symmetries that we will try to implement in the case of gravity 
with anisotropic scaling for general values of z. 



The restoration of general covariance characterized by (L8) still maintains the special 
status of time, keeping it on a different footing from space. We view the fact that "time is 
different" as a virtue of this approach: Indeed, we are looking for possible concessions on the 
side of general relativity that would make it friendlier to the way in which time is treated in 
quantum mechanics, without changing too much of its elegant geometric nature. 

1.2 The minimal theory of gravity with anisotropic scaling 

Here we review the basics of the simplest version of quantum gravity with anisotropic scaling 
[1,2]. 

We assume that the spacetime manifold M is a differentiable manifold equipped with 
the extra structure of a preferred foliation J^ by co dimension-one leaves S of constant time. 
Our starting point is a theory whose only dynamical field is the spatial metric, represented 
in a coordinate system (x = x*,t) on S by components gij{x,t). In a sense, this is the most 
"primitive" implementation of the idea of anisotropic scaling for gravity: The theory exhibits 
no (time-dependent) gauge symmetries, and consequently the spectrum will contain not just 
the tensor polarizations of the gravitons, but also the vector and scalar graviton modes. 

The kinetic term is written using the generalized De Witt metric^ 

The parameter A is left undetermined by the symmetries of the theory, and plays the role of 
an additional coupling constant. 

This "primitive" theory becomes more interesting when we make it gauge invariant under 
foliation-preserving diffeomorphisms Diff(M, J^), whose generators are 

6t = fit), fe^ = r(t,x). (1.10) 



^Throughout this paper, we use the rather loose notation common in high-energy physics, and refer to 
any one-dimensional Abelian symmetry group as (7(1), regardless of whether or not it is actually compact. 
Moreover, we use the same notation also for the infinite-dimensional, gauge version of the (7(1) symmetry. 

SQijki ^^^ i^g viewed as a metric on the space of symmetric 2-tensors. As in [1,2], we will denote its inverse 
by Qijki, to distinguish it from G*-'* with all four indices lowered via gij. 



The minimal multiplet of fields now contains gij plus the lapse function N and the shift vector 
Ni, which transform under ( l.lOl ) as 



Sgij = dii^gjk + OjC'^Oik + C''dk9ij, 

6N., = d,eN, + edjN, + eg,, + JNi + /iVi, (1.11) 

5N = (,^djN + fN + fN. 

The lapse and shift A^ and Ni , known from the ADM parametrization of the spacetime metric 
in general relativity [4] , play the role of gauge fields of the foliation-preserving diffeomorphism 
symmetry. Indeed, ( 1.11[ ) shows that A^ (or, more precisely, log N), and A^* transform as gauge 



fields under the gauge transformations, 

SN ■ • ■■ 

_ = / + ..., 6N^=e + ---. (1.12) 

(Here the ". . ." stand for the standard Lie-derivative terms in ( |1.11| ).) As a result, it is natural 
to assume that A^ and Ni inherit the same dependence on spacetime as the corresponding 
generators ( |1.10| ): While Ni{t, x) is a spacetime field, N{t) is only a function of time, constant 
along the spatial slices E. Making this assumption about the lapse function will lead to the 
minimal theory of gravity with anisotropic scaling. 

The covariantization of the "primitive" theory is accomplished by replacing the spatial 
volume element with its covariant version, 

VO^^N, (1.13) 

and by trading the time derivative of the metric for the extrinsic curvature Kij of the leaves 
of the foliation J^, 

gij -^ 2Kij = — {g,j - ViNj - VjNi) . (1.14) 

In this fashion, we obtain the minimal realization of the idea of anisotropic scaling in 
gravity. This minimal theory is sometimes referred to as "projectable," because the spacetime 
metric assembled from the ingredients gij, Ni and A^ satisfies the axioms of a "projectable 
metric" on (M, J^), as defined in the geometric theory of foliations. The action of the minimal 
theory is 

5 = A f dtd'^x^N {KijK'^ - XK^ - V) , (1.15) 

where K = g^^Kij. The potential term V is an arbitrary Diff(S)-invariant local scalar func- 
tional built out of the spatial metric, its Riemann tensor and the spatial covariant derivatives, 
without the use of time derivatives. 

Around a free-field fixed point with dynamical exponent z, we will measure the scaling 
dimensions of fields in the units of the spatial momentum: [di] = 1. In these units, the 
volume element in the action is of dimension [dtd^x] = —z — D, suggesting the natural 
scaling dimensions for the field multiplet, 

[g,j] = 0, [Ni]=z-l, [N]=0. (1.16) 



This scaling further imphes that [k^] = z — D. 

If we wish for the theory to be power-counting renormahzable, it is natural to start the 
analysis of possible terms appearing in V at short distances. Around a hypothetical Gaussian 
fixed point, power-counting renormalizability in D-|- 1 dimensions requires V to be dominated 
by terms with 2D spatial derivatives, implying in turn that the dynamical critical exponent 
should be equal to z = D. For example, in 3 -|- 1 dimensions, there is a natural potential 

Vuv = w^CijC'^ + ..., (1.17) 

where t/; is a dimensionless coupling, and C*-' = e*'^^Vfc(i?;^ — ■^RS''^) is the Cotton tensor. In 
( |1.17| ), we have indicated only the part of the potential that is dominant in the ultraviolet, 
with ". . ." denoting the relevant terms which contain fewer than six spatial derivatives and 
become important at longer distances. 

Under the influence of the relevant terms, the theory will flow, until it is dominated at 
long distances by the most relevant terms. In this regime, it makes sense to reorganize the 
terms in V by focusing on those most dominant in the infrared: 

Vm = -ij\R-2A) + .... (1.18) 

Here /x and A are dimensionful couplings of dimensions [/i] = z — 1 and [A] = 2, and the ". . ." 
now denote all the terms containing composite operators of higher dimension compared to 
the displayed, most dominant infrared terms. 

It is useful to note that the algebra of gauge symmetries Diff(M, J^), and their action on 
the fields, can be obtained simply by taking a nonrelativistic reduction of the fully relativistic 
spacetime diffeomorphism symmetry and its action on the relativistic metric gfj,^- As we will 



see in Section 3.1, a natural extension of this procedure to subleading terms in 1/c leads to a 



natural geometric interpretation of the extended symmetries that are the focus of this paper. 

1.3 Comments on the nonprojectable case 

The minimal, projectable theory can be rewritten in the Hamiltonian formalism, with the 



Hamiltonian similar to (1.2), 



H= f d^x [NUq + N'Hi) . (1.19) 

Here T-Lq and T-Li are again functions of the spatial metric and its conjugate momenta. In 
fact, T-Li takes the same form as in general relativity, T-C = —2VkT^'^^, and Tio depends on the 
choice of V. The main conceptual difference compared to general relativity stems from the 
fact that because N{t) is independent of the spatial coordinates x*, it only gives rise to the 
integral constraint J d^xT-Lo = 0. Consequently, compared to general relativity, the number 
of first-class constraints and hence gauge symmetries per spacetime point has been reduced 
by one. 

The first, most naive temptation how to eliminate this discrepancy and get closer to 
general relativity is to restore the full dependence of the lapse function on space and time by 



hand. This option, often referred to in the hterature as the "nonprojectable case" [1,2], can 
be viewed at least from two different perspectives, which lead to different results. 

First, one can follow the logic of effective field theory: Having postulated a multiplet of 
spacetime fields 

ffi,(i,x), Ni{t,x), N{t,x), (1.20) 

we postulate a list of global and gauge symmetries, and construct the most general action 
allowed. In the case at hand, the natural gauge symmetries are the foliation-preserving diffeo- 
morphisms Diff(M, J^). While ( 1.19| ) is invariant under Diff(M, J^), it is not the most general 



Hamiltonian compatible with these gauge symmetries. As was pointed out already in [1,2] 
and further elaborated in [8,9], in this effective field theory approach to the nonprojectable 
theory, all terms compatible with the gauge symmetry should be allowed in the Lagrangian. 
Promoting the lapse function to a spacetime field gives a new ingredient for constructing 
gauge-invariant terms in the action, 

ViN/N, (1.21) 

which transforms under Diff(M, J^) as a spatial vector and a time scalar. Once terms with 
this new ingredient are allowed in the action, the Hamiltonian is no longer linear in N, but 
the algebra of constraints is well-behaved [18]. The constraint implied by the variation of A^ 
is now second-class, and the expected number of propagating degrees of freedom is the same 
as in the minimal theory. 

Another possible interpretation of the nonprojectable theory also starts by promoting A^ 
to a spacetime-dependent field N{t, x). Instead of specifying a priori gauge symmetries, how- 



ever, one can postulate that the Hamiltonian take the form ( 1.19 ), linear in A^ [19]. This step 
must be followed by the analysis of the algebra of Hamiltonian constraints, which determines 
a posteriori whether this construction is consistent, and if so, what is the resulting structure 
of the gauge symmetries. Here the difficulty is in closing the constraint algebra [1, 19-21] (see 
also [22]): For general V, the commutator of 'Ho(x) with ^o(y) is a complicated function 
of all variables, and the requirement of closure is difficult to implement. One interesting 
exception has been found in the infrared limit [23] (see also [19,24]): Adding vr = gijir''^ as 
another constraint closes the algebra, turning vr and T-Lq into a pair of second-class constraints. 
This infrared theory can then be interpreted as general relativity whose gauge freedom has 
been partially fixed. This is a very appealing picture, but the problem is that it cannot be 
straightforwardly extended to the full theory beyond the infrared limit. However, as we will 
see below, the possibility of interpreting vr as an additional constraint in the infrared regime 
as suggested in [23] will be echoed in the generally covariant theory which we present in 
Section ^. 

In addition to the two perspectives just reviewed, there is another option how to close the 
constraint algebra in the nonprojectable theory, and interpret it as a topological field theory.^ 
This is possible when the theory satisfies the detailed balance condition [1,2], i.e., when the 



®This option was pointed out by one of us in [20]; see also Section 5.4 of [1]. 



potential V in ( |1.15|) is of the special form 

1 5W5W 

for some action functional W{gij) which depends only on gij and its spatial derivatives. In 
such cases, it is convenient to introduce a system of complex variables, defined as 

1 6W 1 6W 

a'l = ivr^J + ^, a'^ = -ivr*^ + "-f-. (1.23) 

K^ bgij K^ dQij 

Under the Poisson bracket, these variables play essentially the role of a creation and annihi- 
lation pair, their only nonzero bracket being 

K^(x),a^^(y)] = -^ /'f (1.24) 

The Hamiltonian constraints Tii and Tio can be expressed as simple functions of the 
complex variables, 

W = iVj {a'^ - a'^) , no = ya^^g,,Ma''. (1-25) 

The problematic commutator of ?^o(x) and 'Ho(y) is still rather complicated, but it clearly 
vanishes when a^^ or a*-^ vanish. The constraint algebra can thus be closed by declaring Tii, 
together with either a*-' or a*-' to be the primary constraints. This would then guarantee that 
the original Hamiltonian constraints Tii and Hq, as well as all their commutators, are zero on 
the constraint surface. 

This step can be made more precise as follows. Because o*-' and a^^ are complex conjugates 
of each other, it is not possible to declare only (say) a*-^ to be first-class constraints, at least 
not without making a*-^ and a*-^ formally independent. Instead, we accomplish our goal by 



declaring both a*-^ and a*-' as constraints. Because their commutator ( 1.24 ) is nonzero, these 
constraints are second-class and do not imply any additional gauge symmetry. 

However, a pair of second-class constraints can often be interpreted as a first-class con- 
straint, together with a gauge-fixing condition. We can interpret the theory with the second- 
class constraints a^^ and a^^ in this fashion: First, we choose a*-^ — o*-' as the first-class 
constraint. The gauge symmetry generated by this constraint acts on g^j via 

Sgij = Xij{t,x.), (1.26) 

with Xij an arbitrary spacetime-dependent symmetric two-tensor. This is just the topolog- 
ical gauge symmetry as introduced originally by Witten [25,26], here acting on the spatial 
component of the metric. The theory is then fully specified by the choice of a gauge-fixing 
condition for the topological gauge symmetry. In our case, this choice should restore a*-^ and 
a*-' as second-class constraints. Choosing a^^ + a^^ as the gauge fixing condition is certainly a 
consistent possibility; however, a more interesting scenario is available when we Wick-rotate 
the theory to imaginary time, t = —it. This case is of particular interest because topological 



10 



field theories are typically formulated in imaginary time. In this regime, o*-' and o*-' are now 
real, instead of being complex conjugates. We can then select an asymmetric gauge-fixing 
condition, for example 

1 5W 

This equation is a flow equation for gij as a function of the imaginary time r, reminiscent of the 
Ricci flow equation and its cousins. We have thus obtained a topological field theory associated 
with the flow equations on Riemannian manifolds. Indeed, the number of topological gauge 
symmetries ( |1.26 ) is the same as the number of field components of gij: The theory has no 



local propagating degrees of freedom in the bulk. 

The original Diff(M, J^) symmetry can be viewed as a separate gauge symmetry in ad- 
dition to the topological gauge symmetry ( 1.26 ). However, because the action ( p. . 1 1| ) of 



Diff(M, J^) on gij is a special case of (1.26), including Diff(M, J^) explicitly leads to a redun- 
dancy in gauge symmetries, and triggers the appearance of "ghost-for-ghosts" in the BRST 
formalism. In this respect, the structure of the gauge symmetries is very similar to the 
conventional topological field theories of the cohomological type [25, 26] such as topological 
Yang-Mills theory. 

In this paper, we are interested in gravity with bulk propagating degrees of freedom, 
whose spectrum contains the tensor polarizations of the graviton but not the scalar mode. 
Therefore, we do not pursue the nonprojectable theory further, and look for the missing gauge 
invar iance elsewhere. 

2. Global t/(l)s Symmetry in the Minimal Theory at A = 1 



In general relativity, the value A = 1 of the coupling in ( 1.15 ) is selected - and protected from 
renormalization - by the gauge symmetries Diff (M) of the theory. It is perhaps surprising 
that the case of A = 1 plays a special role in the minimal theory with anisotropic scaling as 
well [1,2]. One can see this by examining the spectrum of the linarized fluctuations around 
the flat space solution. 

For simplicity, we will now assume that the flat spacetime 

gij = 5ij, Ni = 0, N = l (2.1) 

is a solution of the equations of motion, and expand the metric to linear order around this 
background, 

gij = Sij + Khij, Ni = KUi, N = 1 + Kn. (2.2) 

Since n is not a spacetime field but only a function of time, its equation of motion gives 
one integral constraint, and does not affect the number of local degrees of freedom or their 
dispersion relations. Therefore, we only consider the equations of motion for the spacetime 
fields hij and rij. 
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It will be convenient to further decompose the hij and rij fluctuations into their irreducible 
components, 

hij = Sij + diWj + djWi + ( didj - -^^ijd'^ \ B + —6ijh, (2.3) 

where the scalar h = ha is the trace part of hij, while Sij is symmetric, traceless and transverse 
(i.e., divergence-free: diSij = 0), and Wi is transverse; and similarly, 

ni = Ui + diC, (2.4) 

with Ui transverse, diUi = 0. It is also useful to decompose the linearized gauge transforma- 
tions, 

r(x,t) = Ci(x,t)+5ir?(x,i). (2.5) 

In this decomposition, (^* satisfy did = 0) ^-^d therefore represent the generators of linearized 
volume-preserving spatial diffeomorphisms. The linearized gauge transformations act on the 
irreducible components of the fields via 

6sij = 0, 6wi = Oi SB = 2r], 5h = 25^r/, 

Sui = Ci, 5C = i,. (2.6) 

These rules suggest a few natural gauge-fixing conditions. For example, we can set Uj = and 
C = 0, which leaves the residual symmetries with time-independent Ci{^) ^-nd r/(x), or Wj = 
and -B = 0, which fixes the gauge symmetries completely. In either gauge, the spectrum of 
linearized fluctuations around the flat background contains transverse traceless polarizations 
Sij which all share the same dispersion relation (dependent on the details of V) , and a scalar 
whose dispersion is dependent on A. In the vicinity of A = 1, the dispersion relation of the 
scalar graviton exhibits a singular behavior, 

a;2 = (A-l)F(k2,A), (2.7) 

where -F(k^, A) is a regular function of A near A = 1, whose details again depend on V. Thus, 
the scalar dispersion relation degenerates to o;^ = in the limit of A — )• 1 [1,2]. 

2.1 Symmetries in the linearized approximation around flat spacetime 

The spectrum of linear excitations around the flat spacetime shows that the relativistic value 
A = 1 is special even in the nonrelativistic theory, as indicated by the dispersion relation 
of the scalar graviton mode ( |2.7D which degenerates as A — >• 1. This singular behavior was 
explained in [1]: At A = 1, the linearized theory with A = 1 enjoys an interesting Abelian 
symmetry, which acts on the flelds of the minimal theory via 

Srii = dia, 6hij =0, 6n = 0. (2.8) 

Here the parameter a(x) is an arbitrary smooth function of the spatial coordinates, constant 
in time: 

d = 0. (2.9) 
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Since the generator a is independent of time, it is natural to interpret this infinite-dimensional 
Abelian symmetry as a global symmetry: In the nonrelativistic setting, it is the hallmark of 
gauge symmetries in the Lagrangian formalism that their generators are arbitrary functions 
of time. In order to indicate that the Abelian symmetries generated by a(x) represent a 
collection of U{1) symmetries parametrized by the spatial slice S of the spacetime foliation, 
we will refer to this infinite-dimensional symmetry by U {!)■£■ 

At this stage, it is interesting to note that the C^(l)s symmetry looks very reminiscent 
of a residual gauge symmetry in a gauge theory, in which some sort of temporal gauge has 
been chosen. As we will see in the rest of the paper, this intuition is essentially correct, but 
the specific realization of this idea in the full nonlinear theory will be surprisingly subtle. 

In order to see that U{1)t, is indeed a symmetry of the linearized theory at A = 1, it is 
instructive to restore temporarily A, and evaluate the variation of the action under (|2.8| ) in the 
linearized approximation (which we denote by "~"), while allowing a to be time dependent: 

5„5« 2 [dtd^'^U f£-l{d^fB + ^^^^d'h] -2a{X-l){dyc\ . (2.10) 

At A = 1, the last term drops out, and the action is invariant under time-independent a. 
Note also that the term proportional to a in ( ^.10 ) is not gauge invariant under (2^), unless 
A = 1 when it equals 

^~'^ {{d^fB-d^h}^R, (2.11) 



D 

which we recognize as the linearized Ricci scalar of gij . 

Given this global f^(l)s symmetry, it is natural to ask whether it can be gauged. At 
A = 1, this process can be easily completed in the linearized theory. We promote a to 
an arbitrary smooth function of x and t, introduce a gauge field A{-x,t), and postulate its 
transformation rules under the gauge transformations, 

6aA = a. (2.12) 

The gauging is accomplished by augmenting the action by a coupling of A to the linearized 
Ricci scalar, 

S^ = -^^^^ I dtd^'^A {{d^fB - d^h} . (2.13) 

It is easy to see that in the gauged theory, the scalar mode of the graviton has been eliminated 
from the spectrum of physical excitations: With A = as our gauge choice, the equations 
of motion are the same as in the original theory with the global U{1)y. symmetry, plus the 
Gauss constraint 

{d'^fB-d'^h = {). (2.14) 

This Gauss constraint eliminates the scalar degree of freedom, leaving only the tensor modes 
of the graviton in the physical spectrum of the linearized theory. 
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2.2 The nonlinear theory 

We would now like to extend the success of the U{1) gauging from the linearized approxima- 
tion to the full nonlinear theory. Before we can proceed with the gauging, however, we must 
first check whether U{l)j] extends to a global symmetry of the nonlinear theory. 

In the linearized theory before gauging, the parameter a(x) of the infinitesimal f7(l)s 
transformation was independent of time, and consequently we interpreted ?7(l)s as a global 
symmetry. In the nonlinear theory, the linearized transformation of rii in ( p.8D simply be- 
comes^ 

S^Ni = NViU. (2.15) 

However, the condition ( |2.9[ ) expressing the time independence of a is not covariant under 
Diff(M, J^). The correct covariant generalization takes the following modified form, 

a - N'ViO = 0. (2.16) 

This condition of vanishing covariant time derivative of a is indeed invariant under Diff(M, J^). 
In the full nonlinear theory, the gauge field will transform as a spatial scalar and a time 
vector under Diff(M, J^), 

6A = fA + fA + Cd^A, (2.17) 

and the gauge transformation of the gauge field becomes 

6aA = a-N'Via. (2.18) 

It follows from (2.18) and ( |1.16| ) that the scaling dimensions of a and A are given by 

[a]=z- 2, [A] = 2z-2. (2.19) 

In the process of evaluating the variation of the action under the general a transformation, 
we will encounter a particular combination of the second spatial derivatives of cjij , which can 
be expressed as the trace of the time derivative of the Ricci tensor: 



g'^R^j = [g^'^g^' - g'^ g^') V^V^gM- (2.20) 

This formula also implies 

(VaR)- = -Va [r'' - ^^5*^) 9^, + Va {g'^g'' - g''g'') v.v.^h- (2.21) 



^The explicit multiplicative factor of A'^ in (2.15) is explained by the requirement of matching the tensorial 
properties of both sides in (2.15) under Diff(M, J^). Thus, it is in fact Ai = Ni/N that transforms as the 
spatial projection of a spacetime vector field under Diff (A/, T) and as the spatial part of a gauge potential 
under U(l), with 5Ai = ViCt. 
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It is now straightforward to see that there is an obstruction against extending the global 
symmetry to the full nonlinear theory, at least in dimensions greater than 2 + 1. Indeed, for 
the variation of the action we get 

^aS = -^ fdt d^x ^ {gij - ViNj - VjNi) [g'^g^' - g'^ g^'^ (VfcV.a + V^V^a) 

j dtd^'^^a [g'^ij - IG'^'^'VkVeViN,'^ , (2.22) 



2 



where in the second line we have integrated by parts twice, dropped the corresponding spatial 
derivative terms, and used ( 2.20| ). The last, triple-derivative term in ( 2.22 ) can be simplified 



usmg 

which yields 

SaS = -^ fdtd^^^a {g'^Ri, - 2Vj{R^''Nk)] . (2.23) 

Finally, after using the contracted Bianchi identity in the second term, integrating by parts 
in both terms, using ( 2.21 ) and dropping the total derivatives, we obtain 

S^S = ^ f dt d^x V5 (a - N'Vio) R 

-^ [dtd^'^^a (r'^ - ^RgA {gij - V,Nj - VjN,) . (2.24) 



The first line in ( p. 24 ) vanishes for the covariantly time- independent a by virtue of ( p. 16 ) 



but the second line represents an obstruction against the invariance of S, even when a is 
restricted to be covariantly time- independent. 

Another way of seeing the origin of the nonlinear obstruction against the ?7(l)s invariance 
of the minimal theory is the following. On the components Ni of the shift vector, the ^(l)i; 
transformations act as gauge transformations on the components of an Abelian connection. 
Define 

Fij = diNj - djN^. (2.25) 

Clearly, this is the field strength of Ni interpreted as a connection associated with {7(l)s- 
Thus, Fij are invariant under U {!)■£, and transform as components of a two- form under 
Diff(S). However, Fij do not transform as two-form components under time-dependent spa- 
tial diffeomorphisms. 

In this new notation, the action with A = 1 can be rewritten as 

5 = ^ / rfi rf'^x ^ {g,, {g^^g^' - g^^ g^') g^ - F'^F,, 

- AR'^NiNj + AN' {V^^j - g'^ViCjjk) ]-^ fdtd^jc^ NV. (2.26) 
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While the first two terms in ( p. 261) and the potential term are manifestly invariant under 
U{1)y:, the terms with explicit factors of A'^j - which are required by the requirement of 



Diff(M, J^) invariance - are not, and their variation reproduces ( |2.24| ). Intuitively, the ob- 
struction can be related to the fact that Ni plays a dual role in the theory. First, as we have 
seen in ( |1.12| ), Ni is the gauge field of the time-dependent spatial difFeomorphisms along S. 
The second role is asked of Ni/N in our attempt to extend the gauge symmetry, and make 
Ni/N transform as the spatial components of a U{\) gauge field. 

3. Gauging the f/(l)s Symmetry: First Examples 

Our intention is to gauge the action of the global U{1)y, in the minimal theory with A = 1. As 
we have seen in Section 4.1.3| , such gauging is possible in the linearized approximation, and 



it has the desired effect of eliminating the scalar polarization of the graviton. However, in 



Section 2.2 we found an obstruction which prevents t/(l)s from being a global symmetry of 
the minimal theory at the nonlinear level, and therefore precludes its straightforward gauging. 
More precisely, we found that the variation ( 2.24| ) of the action under an infinitesimal U{1) 



gauge transformation a(t,x) consists of two parts, 

— I dt d^x V5 (a - N'Vio) R (3.1) 



and 



-^ fdtd^x^a (r'^ - ^RgA {(ji, - ViNj - VjN) . (3.2) 

If C/(1)e were a global symmetry, the first step of the Noether procedure would be to add a 
Noether coupling term to the action, 

Sa = -^ I dtd^i^y/gAR. (3.3) 

The U{1) variation ( |2.18| ) of ^ in ( |3.3D indeed cancels ( |3.1| ). However, at this stage the Noether 
procedure breaks down, and ( |3.2D represents the obstruction against gauge invariance. 

In order to allow a straightforward gauging on f7(l)x;, we will have to find a mechanism 
which eliminates this obstruction. Before proceeding with that, we first examine the geometric 
origin of C/(l)s as a natural gauge symmetry, and consider several illustrative cases in which 
the gauging can be completed because the obstruction automatically vanishes. These include 
a generally covariant nonrelativistic gravity theory in 2 -|- 1 dimensions, and an interacting 
Abelian theory of gravity in general dimensions. 

3.1 Geometric interpretation of the U{1) symmetry 



The transformation rules (1.11) of Diff(M, J^) on the gravity fields can be systematically 
derived [1,2] from the action of relativistic diffeomorphisms Diff (M) on the spacetime metric 
g'^j/, as the leading order in the nonrelativistic 1/c expansion. 
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It was already noted in [2] that the gauge field A{t, x) and the U{1) symmetry - which we 
introduced in a rather ad hoc fashion in Sections |^ and ^ - both acquire a natural geometric 
interpretation in the framework of the 1/c expansion: It turns out that A is simply the 
subleading term in the 1/c expansion of the relativistic lapse function, and ^(1) corresponds 
to the subleading, linearized spacetime-dependent time reparametrization symmetry of the 
relativistic theory. In this section, now make these observations more precise. 

In Section IL^, we reviewed some of the difficulties faced in the attempts to promote the 
lapse function to a spacetime field, 

A^(t) ^iV(t,x). (3.4) 

In physical terms, the attempts to restore A^ as a spacetime field can be motivated by the 
desire to restore the information carried by the Newton potential in general relativity (for 
generic gauge choices). The geometric understanding of the gauge field A and the gauge 
symmetry U{1) shows how the generally covariant theory with U{1) t< Diff(M, J^) symmetry 
restores the Newton potential and avoid the difficulties of the nonprojectable theory: Instead 
of promoting N{t) into a spacetime field as in ( |3.4[ ), we keep N{t) as the leading term of the 
lapse, and introduce the subleading term A{t, x) in the 1/c expansion, at the order in which 
the Newton potential enters in the nonrelativistic approximation to general relativity: 

N{t)^N{t)-\A{t,x). (3.5) 

c^ 

Thus, it is only the subleading part of lapse that becomes a spacetime field. 

It is useful to stress that the 1/c formalism of this section is just a trick, whose sole 
purpose is to provide a geometric explanation of the action of U{1) i< Diff(M, J^), by taking 
the formal c — )• oo limit of the relativistic DifF(M) symmetry. The "speed of light" c is a formal 
expansion parameter, and should not be confused with the physical speed of light which will 
be generated in our theory at large distances as a result of the relevant deformations. 

3.1.1 The gauge field A and the Newton potential 

In order to reproduce the field content of the theory, and the transformation rules under the 
gauge symmetries, we start with a relativistic spacetime metric, and expand it in the powers 
of 1/c as follows: 



r, NiN' 2NA Ni \ 

c^ c^ c 

— + ..., gij + ... / 



(3.6) 



This step is complemented by a similar expansion of the relativistic spacetime diffeomorphisms 
with generators C'^, 

C° = c/(x,t)--^%^ + ..., C = r(x,t) + .... (3.7) 

c iV 
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In both cases, ". . ." refer to terms suppressed by 1/c^ compared to those displayed. In the 
transformation rules, the derivative with respect to the relativistic time coordinate is written 
as d/dx^ = {l/c)d/dt; it is then the nonrelativistic time t which is held fixed as c — )■ cxd. 

Taking the c — )• oo limit first requires dif = 0, which means that the infinitesimal time 
reparametrizations f{t) are restricted to be only functions of time. In addition, in accord 
with (|3.5|), we insist that N be only a function of time. The transformation rules (1.11) under 



the foliation-preserving diffeomorphisms then follow from the c — >■ oo limit of the spacetime 
diffeomorphism symmetry. In addition, we get 

n) = ^^"' 

6aA = a- N'Via, (3.8) 

with all other fields invariant under 5a- We see that the U{1) gauge symmetry of interest is 
geometrically interpreted as the subleading part of time reparametrizations in the nonrela- 
tivistic limit of spacetime diffeomorphisms in general relativity. 

This embedding of the gauge field A into the geometric framework of the 1/c expansion 
sheds additional light on the physical role of A in the theory. Recall that in the leading order of 
the Newtonian approximation to general relativity, the goo component of the spacetime metric 
(in the natural gauge adapted to this approximation) is related to the Newton potential $ 
via 

500 = - (l + ^2$ + . . . y (3.9) 



Comparing this to ( |3.6| ) , we find that our gauge field A effectively plays the role of the Newton 
potential, 

yl = -$ + .... (3.10) 

As we will see in Section (^]^), this relationship is corrected by higher order terms already at 
the next order in the post-Newtonian approximation. 

3.1.2 Extending the 1/c expansion 

We can also keep the subleading terms in the spatial metric, replacing 

g^,^g^J-^^^ + ... (3.11) 

in ( 13.61 ). Following the rules of transformation for the spatial metric to one higher order in 
1/c^ than before, it turns out that Aij also transforms under a, 

6aAij = aijij + NiVjO + NjVia. (3.12) 

This transformation property of Aij is just what is needed to remedy the noninvariance of 
our action under the local U{1) transformations, by introducing a new coupling 

^Jdt d^x ^ Aij (r'^ - ^RgA . (3.13) 
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Interestingly, this term is also "accidentally" invariant under another Abelian gauge sym- 
metry, which acts on the fields via 

6Aij=V^eJ+Vjei. (3.14) 

The variation of all the other fields under the £i symmetry is zero. The total action is invariant 
under ( |3.14| ): The only term in the action which depends on Aij is ( p. 13 ), and its invariance 
under (|3.14|) is a consequence of the Bianchi identity. 



This new gauge symmetry ( tj.14 ) also has a natural geometrical origin: In the process 



^^ = C-^1M2^ + .... (3.15) 



of decomposing the relativistic symmetries in the powers of 1/c, we could have also kept the 
subleading terms in spatial diffeomorphisms, 

1 e*(t,x) 

The c — )• oo limit of the relativistic diffeomorphisms then implies precisely the transformation 
rules (|3l^ ). 

It thus appears that by extending the gravity multiplet to include both the Newton- 
potential A{t,x) and the field Aij{t,x), we succeeded in finding a formulation of gravity 
in which the U{1) ix Diff(M, J^) symmetry of "nonrelativistic general covariance" is realized 
in the full nonlinear theory without obstructions. In addition, we have also seen that this 
extended gravity multiplet has a clear and natural geometric interpretation in the context 
of the 1/c expansion. These features make this extended theory potentially attractive, but 
a closer inspection shows that the number of propagating degrees of freedom has been once 
again reduced to zero - the theory turns out to be effectively topological. Consequently, the 
spectrum of bulk gravitons in the low-energy limit will not match the prediction of low-energy 
general relativity. 

In order to see this, and to count reliably the number of degrees of freedom, we turn 
once more to the Hamiltonian analysis [10]. Because the subleading fields A and Aij that we 
kept in the 1/c expansion appear in the action without time derivatives, they will all lead to 
constraints in the Hamiltonian formulation of the theory. The full phase space is parametrized 
by fields Ni, A, Aij, gij and their canonical momenta P*, Pa, P^-' and vr*-', implying that 

dimP = 2(D + l)2 (3.16) 

per spatial point. ^ The vanishing of the momenta conjugate to A, Aij and Ni represents 
{D+2){D+l)/2 primary constraints. The condition that the primary constraints be preserved 
in time yields secondary constraints: Insisting on Pa = requires the vanishing of R, and 
similarly P*-' = requires the vanishing of W^ — ^Rg^K In addition, as in general relativity, 
P* = requires Tii = 0. 

Naively, there are thus D{D + 3)/2 + 1 secondary constraints Tii, R and R^^ — ^Rg^^ . 
However, these are not all independent: R^^ — oRg^'' satisfies the Bianchi identity, and R is 



*There is also the canonical pair consisting of N{t) and its conjugate momentum Po{t), which only yields 
an integral constraint and can be dropped for the purpose of counting the local degrees of freedom. 
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proportional to the trace of EJ'^ — hRg^^ , leaving D{D + l)/2 independent secondary con- 



2 

straints. 



All the primary and secondary constraints are first-class: Their commutators vanish on 
the constraint surface. As a result, we have the total number 

Ci = (D + 2){D + l)/2 + D{D + l)/2 = (D + 1)^ (3.17) 



of first-class constraints. This implies, together with ( p. 16 ) and invoking ( |1.1| ), that the total 
number of local propagating degrees of freedom is 

AA = -(dimP-2Ci) = 0. (3.18) 

The theory is effectively topological. 

Since our primary inerest in this paper is to find a theory whose spectrum of gravitons 
matches general relativity at long distances, we will not pursue the extended theory in which 
the gravity multiplet contains the Aij fields, and set Aij = from now on. 

3.2 Generally covariant nonrelativistic gravity in 2 + 1 dimensions 

In 2 -|- 1 dimensions, the Einstein tensor Rij — 2^9ij of the spatial metric vanishes identi- 



cally, which means that ( p.2| ) is zero, and there is no obstruction against gauging the U{l)j] 
symmetry in the full nonlinear theory. The Noether procedure terminates after one step and 
leads to the following action, 

S = ^ f dtcfx^lN {KijK'^ - XK'^ - V) - Ar\ . (3.19) 

This action exhibits the [/(I) k Diff(M, J^) gauge symmetry of nonrelativistic general covari- 
ance in 2 + 1 dimensions, for any choice of V. 

The extended gauge symmetry eliminates the scalar degree of freedom of the graviton. 
To see that, it is convenient to select yl = as the gauge choice. In this gauge, the equations 
of motion are the same as in the minimal model with A = 1, with the addition of the Gauss 
constraint 

i? = 0. (3.20) 

It is this additional constraint which eliminates the scalar degree of freedom of the minimal 
theory. Moreover, since in 2 + 1 dimensions the scalar graviton was the only local degree of 
freedom, the generally covariant theory with the extended U{1) k Diff(M, J^) symmetry has 
no local propagating graviton polarizations. In this sense, it is akin to several other, much 
studied models of gravity in 2 + 1 dimensions, such as standard general relativity or chiral 
gravity [27,28]. 

Because of the absence of physical fluctuations, the geometry of classical solutions in 
this theory can be expected to be quite rigid, just as in the case of its relativistic cousins in 
2 + 1 dimensions. In particular, the Gauss constraint ( |3.20| ) forces the two-dimensional spatial 
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slices to be flat. It is natural to look for deformations of this theory which would at least 
replace the Gauss constraint with the more general condition of constant spatial curvature, 
but there appear to be no consistent deformations that could modify the Gauss constraint to 

R = 2n, (3.21) 



with Q a new coupling constant. However, once we learn in Section 4.2 how to gauge the 
C/(l)s symmetry in the general case of D + 1 dimensions, we will also find a mechanism for 
turning on this new coupling 0. 

3.3 Self-interacting Abelian gravity 

Another way to eliminate the obstruction against gauging t/(l)s is to linearize the gauge 
symmetries of the minimal theory. The fields in the theory with linearized gauge symmetries 
are hij, ni and n. The gauge transformations ^j(t,x) and /(t) act via 

5hij = diij + djii, 

5n, = i, (3.22) 

6n = /, 

and represent the Abelian contraction of Diff(M, J^). 
The kinetic term (with A = 1) takes the form 

Sk = :j dt (Pyi (hij - diUj - djuA {Sik6ji - 6ij6ke) (hke - dkUi - deukj . (3.23) 

In this theory, the obstruction (p^) against gauging vanishes identically, as was already es- 



tablished in our analysis of the linearized approximation to the minimal theory in Section 2.1 
At first glance, it would thus seem that keeping only the linearized gauge symmetries 

would reduce the model to the nonintreracting Gaussian theory studied in Section |2.1| , but in 

fact it is not so. Even though the kinetic term takes the Gaussian form ( p. 231 ), the potential 

term need not be Gaussian. 

Suitable terms in V are integrals of local operators, which are either invariant under 

( |3.22| ), or invariant up to a total spatial derivative. The building blocks that can be used 

to construct such operators are the linearized curvature tensor of the spatial metric, and its 

derivatives. We will denote the linearized Riemann tensor by 

Lijki = 2 (djdkhii - djdghik - didkhje + didihjk) , (3.24) 

and similarly the linearized Ricci tensor by Lij = Likjk and the Ricci scalar by L = L^. 
Clearly, there is an infinite hierarchy of sutable operators, which reduces to a finite number if 
we limit the number of spatial derivatives by 2z. For interesting values of z > 1, the general 
V built from such terms will not be purely Gaussian, leading to a self-interacting theory. 

Thus, in the context of gravity with anisotropic scaling, linearizing the gauge symmetries 
does not necessarily make the theory noninteracting - we find a novel interacting theory of 
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Abelian gravity instead. Curiously, a similar phenomenon has been observed in the case of 
general relativity 25 years ago by Wald [29], where it was shown that by taking the action 
to contain higher powers of the linearized curvature, one can construct a self-interacting 
theory of spin-two fields in flat spacetime with linearized spacetime diffeomorphisms as gauge 
symmetries. In the relativistic case studied in [29], this construction leads inevitably to higher 
time derivatives in the action, and therefore problems with ghosts in perturbation theory. In 
contrast, our nonrelativistic models of self-interacting Abelian gravity do not suffer from this 
problem - their self-interaction results from higher than quadratic terms in V, with the kinetic 
term taking the Gaussian form ( ^.23 ). For suitable choices of the couplings in V, the spectrum 



is free of both ghosts and tachyons. 

For an arbitrary V, the gauging of U{l)j:, is now accomplished by adding a new Gaussian 
term to the action, 

-2 dtd^yiAL. (3.25) 



The theory is now gauge invariant under the linearized action of U{1), 

SA = a, 5ni = dia. (3.26) 



Arguments identical to those in Section 2A show that the theory contains only the tensor 
graviton modes, eliminating the scalar. 

At long distances, the dominant terms in V are those with the lowest number of spatial 
derivatives. Since the only suitable operator with just two derivatives is the quadratic part 
of the spatial Einstein-Hilbert term, 

dtd^x UijLij - ^huLJ , (3.27) 

the theory becomes automatically Gaussian at long distances, and approaches a free infrared 
fixed point with z = 1. This behavior can be avoided if we insist that all operators V are 



integrals of gauge-invariant operators: Since (|3.27D is only invariant up to a total derivative, 



it does not belong to this class. In these restricted theories, the infrared behavior will be 
controlled by Gaussian terms with z > 2. In fact, such self-interacting Abelian gravity 
theories, approaching free-field Lifshitz-type fixed points z > 1, have been encountered in the 
infrared regime of a family of condensed matter models on the rigid fee lattice in [30] . 

While such self-interacting Abelian gravity models might be useful for describing new uni- 
versality classes of bose liquids in condensed matter theory, they do not appear phenomeno- 
logically viable as candidates for describing the gravitational phenomena in the observed 
universe. 

4. General Covariance at a Lifshitz Point 

So far, we focused on the special cases in which the obstruction to the gauging of U{l)j] is 
absent. However, none of the resulting theories of gravity with extended gauge symmetry 
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discussed in Section |3| appear phenomenologically interesting as models of gravity in 3 + 1 
dimensions. 

Here we change our perspective, and present a robust mechanism which allows f7(l)s to 
be gauged in the general spacetime dimension D + 1. This will lead to a theory of gravity 
with nonrelativistic general covariance which reproduces many properties of general relativity 
at long distances. 

4.1 Repairing the global C/(l)s symmetry 



In Section 2^, we found an obstruction that prevents the C^(l)s from being a global symmetry 
of the full nonlinear theory for D > 2. Leaving aside the possibility that this obstruction 
could be cancelled by quantum effects (perhaps by a mechanism similar to [31,32]), we look 
for a way to repair the U{l)s symmetry at the classical level. 

4.1.1 The Newton prepotential 

In order to eliminate the obstruction, we introduce an auxiliary scalar field i', which transforms 
under f/(l)s as 

Sal' = a. (4.1) 

We will refer to this field as the "Newton prepotential." The scaling dimension of v is the 
same as the dimension of a, 

[u] = z- 2. (4.2) 

We can now repair the U{1)y; symmetry by adding a new term to the action, 

+ ^ fdtd^x^Nu (r'^ - ^RgA ViVjU. (4.3) 

The variation of i' in the linear term compensates for the noninvariance of the original action 
of the minimal theory. The term quadratic in z^ is in turn required to cancel the variation of 
Ni in the term linear in u. 

4.1.2 Relevant deformations 

We can check by linearizing around the flat background that the number of propagating 
degrees of freedom has not changed by the introduction of the Newton prepotential terms in 
the action. It is rather unconventional that in the expansion around the flat spacetime, the 
new field i^ enters the action at the cubic order in small fields, i.e., its presence does not affect 
the propagator. This issue is eliminated by noticing that a new term, of lower dimension and 
also invariant under the global U{l)s symmetry, can be added to the action: 

Sn = ^ dt d^x ^ ug'^ {g,, - V,Nj - VjN,) + ^ dt d^^^NuAu. (4.4) 
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Here is a coupling constant of dimension [0,] =2. With the addition of this relevant term, 
the Newton prepotential enters the linearized theory, at the quadratic order in fields around 
the flat spacetime. 

The Newton prepotential enters the action with global f7(l)s symmetry quadratically, 
and can be integrated out by solving its equation of motion, 

Q'^V^VJU + &^K^j = 0, (4.5) 

where we have introduced 

Qij ^ pdj _ \jigij + ^gij_ (46) 



Integrating out u would result in a nonlocal action, because ( [4.5D is solved by 

^o(x,t) = - /"d^x'--l— (x,x')(G'=^E:,,)(x'). (4.7) 

Here we have assumed that the operator G^-'VjVj is invertible, and denoted its Green's 
function by (0*-'VjVj)^^(x, x'). We will not try to determine the exact conditions under 
which this assumption is true. However, for example near the flat spacetime geometry (|2.1D, 



we have 0*-^ V^Vj = OA + 0{hij), where A = 9^ is the flat-space Laplacian. This operator is 
invertible at least in perturbation theory, as long as we keep fi nonzero.^ 

Note that the Green's function is still a local function in t. Note also that the expression 
( [4. 7]) for vq has the right form in order for the action to be f/(l)s invariant after u has been 
integrated out. In particular, 

<5„i/o(x,t) = -J /d^x'— -^— (x,x')G'=^(-VfcV,a- V£V,a)(x') 

= y"d^x'^-^(x,x')(e'=VfcV,«)(x') = a(x,t). (4.8) 

The nonlocality of the action obtained by integrating out v is relatively mild: In particular, 
this nonlocality is purely spatial, along the leaves of the spacetime foliation J^. Such non- 
localities are quite common in rather conventional condensed matter systems. Nevertheless, 
in the rest of the paper, we will keep the action manifestly local, by keeping the Newton 
prepotential z^ as an independent fleld instead of integrating it out. 

4.1.3 Linearized theory with global t/(l)E around flat spacetime 

First we will check that our repair of the global f/(l)s symmetry has not changed the count 
of the number of degrees of freedom. Even with the coupling turned on, the flat spacetime 
geometry 

gij = Sij, Ni = 0, N = l, z/ = (4.9) 



^Another interesting example, which will be important below when we gauge the C/(1)e symmetry, is the 
case in which gtj is the metric of a maximally symmetric space satisfying R = 2$!. In this reference background, 
we have 0'-' ViV-, — 2QA/D (with A the Laplace operator of gij), which is also invertible for Q ^ 0. 
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is still a classical solution of the theory (with A = 0), and we can expand around it. 
The V equation of motion is 

292(z^-C) + /i = 0. (4.10) 

The momentum constraints give 

d^iwi - Ui) = (4.11) 

and „ 

D- 1 



2^d'^v+^-—^[(d'^fB-d'^h\=^. (4.12) 



Setting B = ^ and Wj = fixes gauge completely, and implies (with appropriate boundary 
conditions at infinity) that Uj = and 

2Viv = — ^h. (4.13) 

In this gauge, the remaining equations of motion are 

-Sij + ^^^6ijh + 2{didj - 6ijd'^)C - indijO - ^ = O, (4.14) 



where V2 denotes the quadratic part of V in the linearized theory. Using (4.13), we see 
that the v term cancels the h term exactly, allowing one to determine /i as a function of C, 
and substitute back into ( 4.1C| ). The resulting equation determines the dispersion relation of 



the scalar polarization of the graviton. For example, when we set the cosmological constant 
A = 0, the potential term will be dominated at long distances by V = — i?, and we get 

^ = -dh,j + ^^^{d^dj - 5ijd^){d^B - h). (4.15) 

The metric equation of motion then implies that 



2D ■ 
h = -=^—C, (4.16) 

2-D ' ^ ' 



and the v equation of motion gives 



^^-—^d'^C + DC + {D-2)d'^C = Q. (4.17) 

The spectrum thus contains the transverse, traceless polarizations Sjj with dispersion 

6^2 = k^, (4.18) 

plus a scalar graviton (described in this gauge by C) which exhibits the dispersion relation 

implied by ( [4.17| ), 

2 (Z)-2)k2 
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Note that the scalar mode is inevitably tachyonic at low energies. This is implied by 
the choice of sign in V, determined from the requirement that the tensor polarizations have 



the correct-sign dispersion relation ( 4.18 ). This tachyonic nature of the scalar mode is not 
a cause for any concern, because the model discussed here represents only an intermediate 
stage of our construction - our intention is to gauge the U{l)s symmetry, which will turn the 
scalar graviton into a gauge artifact. 

Note also that taking the regulating dimensionful coupling Q to zero reduces the scalar 
dispersion relation ( [4. 191 ) correctly to the singular dispersion w^ = 0, observed at A = 1 in 



the minimal theory in [1] and in (2.7). 



4.2 Gauging the U{l)j] symmetry 

Having repaired the global C/(1)e symmetry, we can now gauge it. The Noether method 
closes after just one step; adding 

SA,n = -^ fdtd^x^A{R-2n) (4.20) 

to the action makes the theory gauge invariant under the U{1) symmetry. This procedure 
results in the following action of the generally covariant theory of gravity with anisotropic 
scaling, 

S = ^ I dt d^x ^In KijK'^ -K'^-V + u e'^ {2Kij + V^V^z^) -A{R-2a)\, 

(4.21) 
with O*-' a short-hand notation for 

Qij ^ j^ij _ \gi3Jl + ^gij_ (4.22) 

Note that in the theory with the Newton prepotential i/, the issue about the possibility 
of adding the spatial cosmological constant O, raised in the generally covariant theory in 2-|- 1 
dimensions at the end of Section ^]^, has been resolved by the introduction of the Newton 
prepotential. 

Note also that in addition to the newly introduced gauge field A, the theory contains a 

composite 

N ■ 
a = v- N'ViV + —VvViV (4.23) 

which also transforms as a gauge field under the U{1) gauge transformations, 

5„o = d - N^Via. (4.24) 

Moreover, both A and the composite gauge field a share the same transformation properties 
under the rest of the gauge group, 

5a = Cditt + fa + fa. (4.25) 
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As a result, the gauged action stays gauge invariant if we replace 

^^(1-7)^ + 70, (4.26) 

with 7 a real coefficient. 

In fact, the composite field a already made its appearance in the theory with the global 

U{1)t: symmetry presented in Section W^. Up to a total derivative, the relevant term ( [4.4D 

can be rewritten as „ 

40 /■ 
Sn = J Vda. (4.27) 



Hence, the substitution ( [4.26| ) in the J dtd 'x.^^A term in ( 4.2C1| ) will just shift the effective 



value of O. Similarly, substituting ( [4.26| ) in the AR term in the Noether coupling ( 4.20 ) 
effectively shifts the coefficients in ( |4.3| ). 

4.2.1 Hamiltonian formulation 

The structure of gauge symmetries can be verified by analyzing the algebra of Hamiltonian 
constraints of the theory. In addition, this analysis will allow us to obtain the precise count 



of the number of propagating degrees of freedom, using formula (1.1). This approach to the 
count of the degrees of freedom is usually more accurate and more reliable than our previous 
analysis of the linearized spectrum around a fixed solution, for at least two reasons. First, 
it is less background-dependent, because it sidesteps the need to linearize the theory around 
a fixed solution. Secondly, because it is valid for the full nonlinear theory, it excludes the 
possible artifacts of the linearized approximation. 

We will set k = 1 to eliminate additional clutter, and denote the canonical momenta 
conjugate to the spatial metric by H*-': 

U'^ = -— = 2 ^5 {K'^ - g'^K + Q'^v) = tt'^ + 2^ Q'^u. (4.28) 

The lower-case vr*-' are reserved for the standard expressions for the canonical momenta in 
general relativity, 

TT^^ = 2^ {K^^ - g'^K) . (4.29) 

The remaining canonical momenta 

PXx,t) = ^, p,(x,t) = §^, P^(x,t)^^, Po(t)^^ (4.30) 

6Ni ov §A 6N 

all vanish, and represent the primary constraints. The Poisson brackets are 



b.,(x,t),H'=^(y,t)] = - Uf<5j + 6m 5{x - y 



1 

2 

[iV,(x,t),P^'(y,t)] = 5^<5(x-y), [Nit),Poit)] = 1, 

[^(x,t),P^(y,t)] =,5(x-y), [i/(x, t),p,(y, t)] = 5(x - y). 



and zero otherwise. 
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In the canonical variables, the Hamiltonian is given by 

1 



H = I d^x <! N 



2^5 



(ff^ - 2^&^u) Gi.M in'' - 2^e'''u 



where 



2NiVjU'^ + 2^A{R-2n)\ , (4.31) 



Gijki = 2 idikgji + gugjk) - jy—rgijgke (4.32) 



is the inverse of the De Witt metric G^^ of (1.9) for A = 1. 



At this stage, the primary constraints are included in the Hamiltonian with the use of 
Lagrange multipliers C/j(x, t), W(x, t), ^/^(x, t), and Uo{t), 

H ^ H = H + f d^x {UiP' + Up^ + UaPa) + UqPq. (4.33) 

The preservation of the primary constraints under the time evolution given by ( [4.33| ) requires 
that the commutators of the primary constraints with H vanish, yielding the following set of 
secondary constraints which are local in space, 

W = [H,P'] = -2S/jW^, (4.34) 

$ = [H,p,] = -4^N&^ViV,u + 4^N&^gijMe'''i^ - 27ve'^g,JM^^^ (4.35) 

^ = [H, Pa] = 2^{R - 2fi), (4.36) 

and one integral constraint 



^M^ 



j d^xUo = [H,Po\ = j d'^xl^ {H'i - 2^Q'iu) Gi^ki (n'^' - 2^Q' 

+ 2^&NiuVjU + 2^v\. (4.37) 

This integral constraint will not affect the number of local degrees of freedom. To avoid 
unnecessary clutter, we concentrate on the analysis of the local constraints, returning to 
( [4.37 ) only at the end of this section. 



Next, we need to ensure that the secondary constraints are preserved in time. The 
momentum constraints "W take formally the same form as in general relativity or in the 
minimal theory of [1,2]. They are indeed preserved in time, albeit in a slightly more intricate 
way than in the minimal theory or in general relativity. In those cases (see the discussion in 
Section 4.4 of [1]), the commutator of T-C with H only gets a contribution from the NkH 
terms in H. The rest of the commutator between the density of the Hamiltonian and W adds 
up to a total derivative, as a consequence of the transformation properties of a scalar density 
under spatial diffeomorphisms. Here, the argument is more subtle, and the commutator 
contains additional terms, 

[H, W] = -VkiN^Ui) - {ViN'')nk - (VV)$ - {V'A)^. (4.38) 
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However, this expression vanishes on the constraint surface, and no tertiary constraints are 
produced at this stage. 

The time preservation of the secondary constraint ^ requires the vanishing of 

[H, $] = 4^N&^ (ViVj - gijMQ''') l(+[H,'^]-Uo^ = 0. (4.39) 

Unhke the conditions for the time preservation of the primary constraints or the Ti^, condition 
( 4.39| ) depends exphcitly on one of the Lagrange multiphers, U. Therefore, setting [H, ^] = 



yields an equation for U, instead of producing an additional, tertiary constraint. Also, because 
the commutator 

[p,(x), cD(y)] = 4^NQ'^ (ViVj - GijkiQ''') <5(x - y) (4.40) 

does not vanish on the constraint surface, pu and $ represent a pair of second-class constraints. 
It now remains to check the condition for the preservation of ^ in time. After a lengthy 
calculation, we get 

[H^ ^] = +NViW - $ - Vi {N'^i) . (4.41) 

This expression vanishes on the constraint surface. Again, no tertiary constraint is produced, 
and the process of generating the full list of constraints stops here. 

One might be tempted to expect that ^ is a first-class constraint, but that expectation 
is false: The commutator of ^(x) and <&(y) does not vanish on the constraint surface. Con- 
sequently, the first-class and second-class constraints are still entangled, and ^(x) itself is a 
mixture of constraints of both classes. In order to disentangle the constraints, we must first 
evaluate 

,5{V5(^-2f])(x)} /,,^ .^^fc, 



[^(x),$(y)] = _4:iivin__^wi ^NG,,uQ''] (y) 



= A^NQ'i [QijuQ''' - ViV,- j 5(x - y). (4.42) 

This is equal, up to a sign, to the commutator of py and <I> which we obtained in ( [4.401) . 
Hence, it is natural to define 

nA = ^+Pv (4.43) 

TIa commutes both with $ and with p^^ and represents a first-class constraint. 

Having identified %a as the final first-class constraint, we can check that it generates the 
correct U{\) gauge transformations on the fields. In the Hamiltonian formalism, the gauge 
transformation generated by a first-class constraint on an arbitrary phase-space variable (j) is 
given by the commutator of (p with the corresponding constraint [10], for example 

<5,,^(x, t) = -[j d^y a(y , t)nA , </'(x, t)] . (4.44) 

One can indeed use this Hamiltonian formula to check that the gauge symmetries implied 
by the first-class constraints reproduce those that we found in the Lagrangian formulation 
above. 
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Given our analysis of the constraints, we can now evaluate the number of degrees of 
freedom. Altogether, the theory has dimP = D^ + 3D + 4 canonical variables per spacetime 
point. These variables are constrained by Ci = 2D + 2 first-class constraints (P*, Pa, W and 
T-La), and C2 = 2 second-class constraints [p^, and <1>). The number of degrees of freedom N 
per spacetime point is then given by formula ( |1 . 1| ) , 

^=\ (dimP - 2Ci - C2) = ^(D + 1){D - 2). (4.45) 

This correctly reproduces the number of tensor [i.e. transverse, traceless) polarizations of the 
graviton va. D + \ spacetime dimensions. 



Returning to the integral constraint (4.37), we note that its commutation relations with 



the rest of the constraint algebra can be read off from the commutators of H obtained above. 
This follows from the fact that, as in general relativity, the Hamiltonian can be written as a 
sum of constraints, 

H = N I d^xn^ + [ d^x {N'Tii + A^) . (4.46) 



Actually, the role of the integral constraint ( |4.37 ) deserves to be investigated further 



It is plausible that in the theory with nonrelativistic general covariance, where the U{1) 
gauge symmetry mimics the role of relativistic time reparametrizations, one can choose not 
to impose the integral constraint on physical states. This would be eqiuvalent to the omission 
of nonrelativistic time reparametrizations 6t = f{t) from the gauge symmetries, effectively 
setting N(t) = 1. If consistent, this construction would lead to a theory of gravity with 
nonzero energy levels even in spacetimes with compact spatial slices S. In fact, this situation 
was already encountered on flat noncompact E in the context of Abelian gravity in [30] . On 
noncompact E, the possibility of relaxing the integral Hamiltonian constraint will be closely 
tied to the structure of consistent boundary conditions at infinity in gravity with anisotropic 
scaling, whose study has been initiated in [33]. 

4.2.2 Linearization around detailed balance 

In principle, our result ( ^.45)) for the number of degrees of freedom A^ can be checked by 
linearizing the theory around a chosen solution, and explicitly counting the number of prop- 
agating polarizations. However, in order to investigate the spectrum of the linearized theory 
after gauging, we cannot use the flat spacetime as a reference background, because it no 
longer solves the equations of motion if 0, is not zero. 

This makes the analysis of the linearized approximation for general values of the couplings 
algebraically tedious, and we will not present it here in full generality. Instead, we content 
ourselves with testing ( 4.45| ) in the simpler case when the theory satisfies the detailed balance 



condition. Hence, we assume that the potential takes the special form 



,, 1 6W 6W 
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and for coiicreteness we choose 

W 






(4.48) 



Before the C/(l)s is gauged, the theory in detailed balance admits a particularly simple static 
ground-state solution, 



9ij = 9ij{y^), 



N = l, 



Ni = 0, 



z/ = 0, 



(4.49) 



where gij is the maximally symmetric spatial metric which solves the equations of motion of 

W, 



Rij — IjRgij + ^W9ij — 0- 



(4.50) 



In order for this background to be a solution of the theory with the extended U{\) x Diff (M, T) 
gauge symmetry, we must set the spatial cosmological constant J7 equal to 



n 



D 



D 



-A 



w 



(4.51) 



For O > 0, the ground-state geometry is the Einstein static universe, with spatial slices 
S = S^ . Conversely, when fi < 0, the ground state is the hyperbolic version of the Einstein 
static universe, with noncompact S. Its curvature tensor satisfies Rij = -jj-Qij and R = 20. 
We now determine the spectrum of linearized perturbations around this class of ground 



state solutions. The analysis closely parallels that of Sections p.l| and 4.1.3, and we will be 
brief. We expand the metric, gij = gij + Khij , and decompose the linearized fluctuations as 



h 



rii 



in {^) and (^: 

Sij + ViWj + VjWi + i ViVj - j^gijA j B + - 

Ui + VjC, 

with V j the covariant derivative of cjij . The u equation of motion is 

and the momentum constraints give 

/-. 20\ 

0, 

0. 



(4.52) 



(4.53) 



A + — ] [Wi- Ui) 



/2Q . D-l,. ;, 4f) 

'^i{^B + ^^{AB-h) + —u 



(4.54) 



To fix the Diff (M, J^) symmetries, we set Wi = B = n = 0. In this gauge, the momentum 
constraints reduce to 



(A + 2n/D)ui = 0, 



Vi 



Aniy-{D-l)h 



0, 



(4.55) 
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which imphes, with suitable boundary conditions, that Ui is not propagating, and that 

D-1 



nu 



-h. 



Plugging this back into ( 4.53 ) yields 

D- 1 



Ah + -h 

4:0, 2 



AC = 0. 



(4.56) 



(4.57) 



Finally, there is the constraint R — 20 = 0, which plays the role of the Gauss constraint 
in our gauge A = 0. Its linearization around our detailed balance background 



fl-ii=»-i^ 



{D - l)Ah + 2nh 







(4.58) 



shows that h is not propagating. Combining ( 4.58| ) with ( 4.57 ) then implies that AC = 0. 
Hence, the only propagating modes are the transverse traceless polarizations of the graviton 
Sij. In particular, the scalar graviton has been eliminated, and the number of physical degrees 
of freedom agrees with the result of our Hamiltonian analysis ( 4.45| ). 



5. Conclusions 

In this paper, we have found a formulation of the theory of gravity with anisotropic scaling 
in which the gauge symmetry of foliation-preserving diffeomorphisms Diff (M, J^) is enhanced 
to the symmetry of "nonrelativistic general covariance," U{1) x Diff(M, J^). 

The advantage of this construction is that it relies only on the structure of the kinetic 
term in the action ( p.. 151 ) (and, in fact, forces it to take the general-relativistic form with 
A = 1), while the form of the potential term V is left unconstrained. Therefore, we can 
consider the scenario proposed originally in [1,2], in which the theory is defined at short 
distances by a z > 1 fixed point (with V dominated by higher-derivative terms), and is then 
expected to fiow under the influence of relevant terms to z = 1 and isotropic scaling in the 
infrared. This classical scenario will of course receive quantum corrections, which could drive 



the theory outside the range of validity of the covariant action (4.21). In the rest of the paper, 
we limit our attention to the possibility that the long-distance physics is still described by 
the same action ( [4.21| ), with V dominated by the most relevant terms (1.18). 

The flrst good news is that, as a result of the extended gauge symmetry, the spectrum 
contains just the transverse-traceless (tensor) modes of the graviton. The scalar graviton 
mode of the minimal theory has been eliminated. In 3 + 1 spacetime dimensions, the elimi- 
nation of the scalar mode has an interesting consequence in the short-distance regime of the 
theory. Recall that in the minimal theory with the potential dominated at short distances 
by the z = 3 term (LIT]), the scalar mode is the sole physical mode that does not get a 



contribution to its dispersion relation from ( 1.17 ), suggesting that terms with z > 3 would be 
required to achieve a UV completion [2]. In the theory with the extended gauge symmetry. 
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the scalar mode is a gauge artifact, all physical modes aquire a z = 3 dispersion relation at 



short distances from (1.17), and no terms with z > 3 are needed. 

The extended gauge symmetry of the theory with nonrelativistic general covariance has 
even more interesting consequences at long distances, because it improves the chances that 
the behavior of our theory can resemble general relativity in this observationally relevant 
regime. We conclude this paper by previewing how our generally covariant theory compares 
to general relativity at long distances, focusing on the case of 3 + 1 spacetime dimensions. 

Note first that even before we take the long-distance limit, the elimination of the scalar 
mode of the graviton is certainly a good sign for the possible matching against general relativ- 
ity at long distances, and so is the fact that the coupling constant A in the kinetic term is now 
frozen by the symmetries of the generally covariant theory to take the relativistic value A = 1. 
As a result, the number and the tensor structure of the gravitational wave polarizations is 
the same as in general relativity. 

In the infrared limit of our theory, the potential V is dominated by the scalar curvature 



and the cosmological constant term, (1.18|). In this regime, the natural scaling is isotropic. 



with dynamical exponent z = 1. The low-energy physics is best represented in rescaled 
coordinates {x^,x^) and in terms of rescaled fields. First, the new time coordinate 

x° = fit (5.1) 

is defined by absorbing the effective speed of light n into the definition of time. Because 
[fj] = z — 1, this implies that [x^] = — 1 = [x*], in accord with the z = 1 scaling. The rescaled 
fields are defined by 

Nj"" = -N,, A'"" = \a. (5.2) 

This rescaling ensures (i) that N^ carries the canonical dimension implied by z = 1, and (ii) 
that the UiX) gauge transformations are given by the standard relativistic formula 

<5A^/^ = dioF^, bJ^^ = doa^^, (5.3) 

with a = a//i. In the rest of the paper, we will drop the "IR" superscripts, and refer to 
the rescaled fields (5.2) in the infrared simply as Ni and A. 



The action of the infrared theory in the infrared variables is 

Sir = j^^ j dx° d^^^[N [K.jK'i - K^ + R-2K) - A{R-2n)^ + . . . , (5.4) 

where ". . ." denotes corrections due to higher dimension operators, as well as the z^-dependent 
terms in ( [4.21 ) which are unimportant for our arguments below. In (5.4), Kij refers to the 



extrinsic curvature tensor in the infrared coordinates, of canonical scaling dimension equal to 
one; and the Newton constant is given by 
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In the remainder of this section, we comment on three issues: The structure of compact- 
object solutions (which will be relevant for solar system tests), the issue of Lorentz symmetry, 
and the nature of cosmological solutions in the infrared regime of our theory as described by 

(r 



5.1 Static compact-object solutions 

To prepare the ground for solar system tests, consider the infrared limit (|5.4| ) and set the 
cosmological constant to zero. Interestingly, as the Schwarzschild black hole turns out to be 
a solution of this infrared theory. In terms of our fields, this solution will be represented by 



2M\ 



-1 



Qijdx'dx^ = M - ^^ J dr^ + r^dnl (5.6) 

/ 2MV^^ 
^ = 1-11 J , iV = 1, Ni = 0, z/ = 0. (5.7) 

It is straightforward to see that this geometry satisfies the equations of motion of our theory 
for f] = 0, which is the appropriate choice if we are interested in asymptotically flat solutions. 
First, the equations of motion contain the condition R = 2il.. With il = 0, this equation 
is indeed satisfied by the spatial slices (^]^) of the relativistic Schwarzschild metric in the 
Schwarzschild coordinate system. The v and Ni equations of motion are also satisfied, because 
the extrinsic curvature Kij vanishes for static backgrounds. 

Finally, to show that the gij equation of motion are also satisfied, we use a simple but 
intriguing argument. Since the same argument generalizes in a useful way to the case of 
nonzero fi = A, and also of arbitrary dimension, we present this more general case. Start 
with static solutions with Kij = 0, and observe that the equations of motion for g^j, N and 
A are identical to the equations that follow from the following reduced action, 



f d^x^{N -A){R-2n). 



(5.8) 



Similarly, for static solutions with Kij = of general relativity in D + 1 dimensions, the 
corresponding equations of motion are those of the reduced Einstein-Hilbert action, 

f d'^x^dl{R-2A), (5.9) 

where 91 is the general-relativistic lapse function. Consequently, if we identify the N and A 
fields with the lapse function ^ of general relativity, 

m = N-A, (5.10) 

we see that static solutions of general relativity are also solutions of our theory in the infrared 
limit. In retrospect, this mapping also explains the form of A in our representation of the 



Schwarzschild metric (5.7). 
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Note that the relationship ( 5.10 ) between the general-relativistic lapse function 9T and 
the N and A variables of our theory reproduces exactly what we would have expected from 
the geometric interpretation of A as the subleading term in the expansion of the relativistic 



goo in powers of 1/c as obtained in ( |3.6D . Indeed, we start by expanding 

900 = -^^ = -{N - Af ^ -N^ + 2NA + ... (5.11) 



Recall now that A in ( 5.11| ) is the infrared rescaled field ( |5.2D , related to the microscopic 
gauge field by a rescaling factor 1/;U^. Using the fact that /x plays the role of the speed of 



light (as we have seen in (^.l])), the two leading terms in ( [5.11 ) match exactly the leading 
two terms in the expansion 



These arguments prove that the Schwarzschild geometry in the Schwarzschild coordinates, 
with the indentification implied by (5.10), is a solution of the infrared limit of our theory. 



with Q = A = 0. However, in the parametrized post-Newtonian (PPN) formalism [34, 35] 
which is typically used in gravitational phenomenology, the compact-object solution is usually 
represented in the isotropic coordinates. In the case of general relativity, this is just a gauge 
choice, a fact which does not extend automatically to alternative approaches to gravity such 
as ours. Showing that the Schwarzschild geometry in the Schwarzschild coordinates is a 
solution of our theory does not imply that it will be a solution when represented in another 
coordinate system, because only those coordinate changes that belong to the gauge symmetry 
of our model will map a solution to a solution. However, because the transformation from 
the Schwarzschild coordinates to the isotropic ones only changes the radial coordinate, 

r = ,(l + ^)', (5.12) 

while keeping t, 9, (p intact, it is a foliation-preserving diffeomorphism, a symmetry of the 
theory. Consequently, the Schwarzschild solution in the isotropic coordinates, represented by 

gijdx'dx^ = M + — j {dp^ + p^dnf) , (5.13) 

A=il + — ] — , iV = l, Ni = 0, u = 0, (5.14) 

is also a solution of the infrared limit of our theory. Expanding this solution to the required 
order in the powers of M/p strongly suggests that in the infrared regime, the f3 and 7 parame- 
ters of the PPN formalism [34, 35] will take the same values as in general relativity, /3 = 7 = 1 . 
This feature is favorable for the solar-system tests of the theory. 

5.2 Lorentz symmetry 

Perhaps the leading challenge in any attempt to make theories of gravity with anisotropic scal- 
ing phenomenologically viable in 3 + 1 dimensions is the issue of restoring Lorentz symmetry, 
at least at the intermediate energies and distances where it has been so well tested experi- 
mentally. In particular, we need a mechanism ensuring that in the corresponding regime, all 
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species of matter (including the gravitons) percieve the same hghtcones and the same effective 
speed of light. In the minimal theory with anisotropic scaling, this issue arises already for 
pure gravity: At generic values of the couplings, the speeds of the tensor and scalar graviton 
polarizations are not related by any symmetry, and are generally different from each other 
already in the short-distance regime. In contrast, our generally covariant theory has only the 
tensor graviton polarizations, all sharing the same speed at all energies; however, the issue 
reemerges when pure gravity is coupled to non-gravitational matter. If the present theory is to 
be phenomenologically viable, its coupling to matter will have to be analyzed in detail. This 
analysis is beyond the scope of the present paper; we only limit ourselves to one observation, 
which may be useful for the future analysis. 

In general relativity, Lorentz symmetry is a global symmetry associated with the isome- 
trics of the Minkowski spacetime. In gravity with anisotropic scaling, we can adjust the 
couplings such that the flat spacetime geometry continues to be a solution. The global sym- 
metries of this solution will then depend on the precise model of gravity with anisotropic 
scaling. 

First consider the case of the minimal theory reviewed in Section (1.2), with the cosmo- 



logical constant tuned to zero. The flat spacetime (2.1) is a solution, but it does not exhibit 



the full global Lorentz symmetry - the Lorentz boosts, generated by 

6t = bix\ 5x' = hit (5.15) 

(with bi a constant vector), are not foliation-preserving diffeomorphisms. In this theory, the 
nonrelativistic analogs of the Killing symmetries of the flat spacetime solution correspond to 
spacetime translations and space rotations - the solution breaks all possible boost symmetries 
spontaneously, and defines a preferred rest frame. 

In contrast, in our generally covariant theory, we can interpret the Lorentz transformation 
( |5.15| ) as a generator of a transformation belonging to the extended symmetry group U{1) x 
Diff(M, J^). More precisely, the Lorentz transformation ( 5.15| ) should be interpreted as a 



composition of an infinitesimal foliation-preserving diffeomorphism and an infinitesimal U{1) 
transformation. Indeed, restoring the factors of c shows that the variation of t in ( ^.151) 
is suppressed by a factor of 1/c^ compared to the variation of x*, and should therefore be 
interpreted as an infinitesimal ^'^(l) transformation with a = bix"^, accompanied in ( 5.15D by 
the infinitesimal foliation-preserving diffeomophism 

St = 0, 6x' = bit. (5.16) 



When interpreted in this way, the Lorentz transformation ( |5.15| ) is a symmetry of the flat 
spacetime geometry represented in our variables by gij = 6ij, N = 1 and A = 0. This does 
not yet imply that all preferred-frame effects are absent in this background: In particular, 
the Newton prepotential i' is not invariant under the Lorentz boosts, and deflnes a preferred 
frame for the flat spacetime, in which 1^ = 0. The flat background is Lorentz invariant only 
to the extent that the effects of the Newton prepotential can be ignored. 
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5.3 Cosmological solutions 

Moving beyond asymptotically flat spacetimes, it is natural to ask whether our theory has 
interesting cosmological solutions. One can start with a given spacetime geometry in general 
relativity, and investigate whether it satisfies the equations of motion of our theory. The 
answer to this question will again depend on the choice of spacetime foliation. 

For example, arguments identical to those used above for the Schwarzschild metric show 
that the static patch of the de Sitter (or anti-de Sitter) spacetime, represented in our variables 

by 

/ A 2\-l 

gijdx'dx^ = il-—-] dr^ + r^dnl (5.17) 

A=l - [l-—-] , iV = 1, Ni = 0, u = 0, (5.18) 



is a solution of our theory if we set fi = A. 

It is encouraging to see that at least in the time-independent foliations, the de Sitter and 
anti-de Sitter spacetimes are solutions of our theory. In standard cosmological applications, 
however, the cosmological principle selects another natural foliation of spacetime, with homo- 
geneous spatial slices and a time-dependent scale factor a{t). On the face of it, it may appear 
difficult to obtain cosmological solutions of our theory with maximally symmetric and time- 
dependent spatial slices: The equation of motion for A plays the role of a Gauss constraint, 
and implies R = 2Q in the vacuum. Assuming the standard FRW Ansatz 

gij = a\t)^ij, N = l, A = 0, Ni = 0, u = (5.19) 

where 'jij is a time-independent maximally symmetric spatial metric, the scalar curvature of 
gij has to be constant in time. Consequently, if the scalar curvature of 'jij is nonzero, the 
cosmological scale factor must be independent of time. 

Of course, if our spatial slices are flat, the Gauss constraint no longer restricts the time 
dependence of the cosmological scale factor. This requires il = 0. The rest of the equations 
of motion will be satisfled by the de Sitter spacetime in the inflationary coordinates, which 
in the FRW Ansatz ( |5.19| ) corresponds to 

a(t) = e^*, lij = S^j. (5.20) 

The reason for this is again simple but illuminating: With fi = 0, the i' equation of motion 
is satisfied when the metric is flat. With u and A both zero, the remaining equations are 
implied by Einstein's equations if we simply identify the relativistic lapse function with our 
N{t), the relativistic cosmological constant with our A, and H with the Hubble constant. 

Thus, we see that the same de Sitter spacetime in two different foliations is a solution of 
the infrared theory for two different choices of the coupling constant, one with fi = A and the 
other with = and nonzero A. Mapping out the general behavior of cosmological solutions 
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as the coupling constants Q and A are independently varied is one of questions left for future 
work. 

In addition, there are at least two ways out of the potential difficulty with solving the 
Gauss constraint for cosmologically evolving spacetimes with maximally symmetric spatial 
slices of nonzero curvature. First, the equations of motion will change in the presence of 
matter. In the full system of equations for gravity and matter, the Gauss constraint is 
expected to be modified by a matter source, whose time dependence can then drive the time 
dependence of the scale factor in the spatial metric. The second possibility is related to the 
gauge freedom we have in describing cosmological solutions in general relativity: Instead of 



the standard FRW ansatz which leads to (5.19), one can choose coordinates in which the 



spatial metric is not only maximally symmetric but also constant in time. When we express 
the FRW geometry in such coordinates, the time-dependent scale factor appears in the dt"^ 
term in the metric, and non-zero components of the shift vector Ni are typically generated. 
In general relativity, this coordinate representation of FRW cosmologies is a legitimate albeit 
slightly unconventional gauge choice. In our theory, this parametrization of FRW universes 
has the advantage of being compatible with the vacuum Gauss constraint R = 20. 
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